Field theory of statics and dynamics of glasses: rare events and barrier distributions 
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We study thermally activated dynamics using functional renormalization within the field theory 
of randomly pinned elastic systems, a prototype for glasses. It appears through an essentially non- 
perturbative boundary layer in the running effective action, for which we find a consistent scaling 
ansatz to all orders. We find that in the statics the boundary layer describes the physics of rare 
low energy metastable states as suggested by a phenomenological droplet picture. The extension 
to equilibrium dynamics mandates a non-trivial broad distribution of relaxation times for the large 
scale modes, hitherto neglected in analytic calculations. 



Extremely slow dynamics and dominance of rare fluc- 
tuations are ubiquitous properties of complex and disor- 
dered materials. Such behavior occurs in a panoply of 
systems involving many interacting degrees of freedom, 
from spin glasses to supercooled liquids to the pinned 
elastic medium. We study the latter here as a prototype 
due to its relative simplicity. It describes numerous ex- 
perimental systems such as interfaces in magnets Jj], ||, 
charge density waves and Wigner crystals [^|, or vortex 
lattices in superconductors [0 . The elastic medium is de- 
scribed by a displacement field $(r), with Hamiltonian 



H[u] 



d d r 



|V$| 2 + V(<I>(r),r) 



(1) 



where V($, r) is a zero-mean Gaussian random potential 
with variance V($, r)V(&, r') = R($-$')5 {d) (r-r') (an 
overline indicates disorder averaging). Different physical 
cases are described by varying longitudinal dimensional- 
ity (d), transverse dimensionality (N — the number of 
vector components of $ - we take N = 1 for simplicity 
here), and The Langevin equation of motion is 



ridt^rt = V 2 $ rt + r) + C(r, t), 



(2) 



with friction r), thermal noise £(r, t) and random force 
/($,r) = -dV(<f>,r)/d$ (denoting $(r,t) = $ rt ). 

A minimal physical description of disordered glasses 
is provided by the droplet scenario j^]. In its simplest 
form, it supposes the existence, at each length scale L, 
of a small number of excitations of size 5$ ~ IS above a 
ground state, drawn from an energy distribution of width 
8E ~ L e with constant weight near SE = 0. Although 
the applicability of the droplet picture to more complex 
systems such as spin glasses is controversial, it appears 
to describe simpler models such as Eq. (|l|) relatively well, 
at least in low dimensions,^ and we summarize its main 
conclusions. Static thermal fluctuations at a given scale 
are dominated by rare samples/regions with two nearly 
degenerate minima. In particular, the (2n) th moment of 
$ fluctuations is expected to behave as 



The droplet picture supposes the long-time equilibrium 
dynamics is dominated by thermal activation between 
these quasi-degenerate minima controlled by barriers of 
typical scale Ul ~ . Little is known about the dis- 
tribution of these barriers, but there is some evidence 
M that ip « 6. Even a modest distribution of barriers, 
however, due to the Arrhenius law tl ~ e UL ^ T , yields 
relaxation time scales with an extremely broad distribu- 
tion as T — > 0. An important unanswered question is the 
behavior of relaxation time moments, defined e.g. by 



dU L p(U L )e 



Ml IT 



(4) 



and the response function R, 
The moment t£ may behave as t£ 



q,t-t' 



with q = L 

dt> q (t)/dC q (t>). 

e a(n)u L it w ^h a(n) > ?i, or grow even faster with 
lnr™ 3> L e . A theory of these timescales is crucial to un- 
derstanding both equilibrium response and correlations 
and to near-equilibrium phenomena such as creep j^, 0]. 

Despite the wide applicability and extensive theoreti- 
cal studies of the elastic model, Eqs. ([jJH), analytic treat- 
ment for d > 1 of the non-zero temperature glassy behav- 
ior has been limited to mean-field (N — > oo) || ^ and 
infinite dimension (d — * oo) limits [fl0[ in which ther- 
mal activation over divergent barriers (Ul ~ L e ) is not 
included. The existence of these growing barriers ap- 
pears to be captured by the field-theoretic Functional 
Renormalization Group (FRG) method 11 [fI (fo r 



8\ T 2nQ 



(3) 



e = 4— d <S^1), extended to non zero temperature |14L fL5l, 
but until now neither the rare events nor fluctuating bar- 
riers have been obtained in this approach. 

In this paper, we investigate how properties reminis- 
cent of droplets can emerge naturally from the FRG. 
Although temperature is formally irrelevant (in the 
RG sense), thermal effects give rise to an intricate 
non-perturbative boundary layer structure. Neverthe- 
less, a consistent scaling of this boundary layer can be 
understood, providing a field-theoretic framework for 
thermally-activated glassy dynamics. This structure ap- 
pears to support solutions of the droplet type, with in 
particular rare thermal fluctuations obeying Eq. (|3|) and 



2 



a distribution of timescales t" described by a non-trivial 
function a(n). Details can be found in 

We begin by considering the statics. The theory is 
formulated in terms of the disorder-averaged partition 
function for n replicas $ = ($1, • • • , 




(5) 



A more general distribution for V(Q, r) (still 
independently distributed at each r) is ex- 
pressed through the characteristic function 
U[$] = Er=2TwE 01 ... tt ^W[$ 01 ,---,$ 0fe ]. Here 
is the fc th connected cumulant of 
F($,r) at fixed r (note ($,$') = R(<f> - $'))■ 



By translational invariance, these cumulants satisfy 
($! + A, • • • , $ fc + A) = S« ($1, • • • , $fe), so t/[$] isa 
fully symmetric, translationally-invariant function of $. 
As is well-known |ll|, zero-temperature power-counting 
shows that all for k > 2 are irrelevant, and the 
second cumulant i?($ — <£>') is marginally relevant just 
below d = 4. One therefore finds that for e = 4 — <i <C 1, 
the fixed point function i?*(<£>) ~ e <C 1 and ( also of 
0(e). This smallness justifies the truncation « for 
fc > 2 at one-loop, to leading order in e at T = 0. 

Prior investigations of thermal effects have presumed 
this truncation remains valid at T > |l4|, 15 . We 
check this assumption by considering the corrections due 
to S^ 3 ' ^ 0. A straightforward Wilson scheme calculation 
keeping only k = 2, 3 terms to one loop gives 



J 



diR{<P) = (e - 4C + ($d<s,)R{$) + T ;j R"($) + S^(0, 0, $) + -R"($) 2 - R"(0)R"($) 



123) 



(26 - 2 - 6C + C^^)^ 3 )($ 123 ) + sym (j^S^ ($ 123 ) + ^R"($ 13 )R"(<i>23)) 



+3sym ( R"($i 2 )(^($ii 3 ) - §{? ($i23)) + R"($i 2 )R"($i 3 ) 2 - ^R"($i 2 )R"($ 23 )R"($3i) 



5(3) , 



(6) 
(7) 



where A; = Ae 1 is the running cutoff, <5>y = $i — $j, 
= <!•,.<!>.,.«!>;.. With A d = l/(2 d - 1 Tr d / 2 T(d/2)) 1 
the rescaled functions ft($) = A ci Af _4+4C J R($A^ c ) and 
^ (3) ($i 23 ) = A 2 A 2d - 6+6 ^( 3 )($ 12 3Ar C ). xhe ra ^ of 
thermal energy T to pinning energy L e defines T; = 
A d TAf, where 6» = d- 2 + 2C>0 indicates the dom- 
inance of energy over entropy as appropriate for a glass. 
It was convenient to define R($) = R(&) - J R"(0)$ 2 /2, 
since i?"(0) does not feed into higher cumulants. 

The naive T > analysis assumes <C eR to con- 
vert Eq. (||) into a decoupled equation for -/?($). Its solu- 
tion yields two regimes. For $ ~ 0(1) ^> T/, the solution 
i?/($) converges to a zero temperature fixed point func- 
tion pH , with a small set of universality classes. Q. [iri|] 
For, e.g. V r ($, a;) periodic on the interval < $ < 1, the 
force-force correlator A($) = -R"($) = Min„ e2 |[(<i> - 
n — 1/2) 2 — 1/12]. Non-analytic behavior at small <E>, 
A T=0 ($) - A(0) ~ -x|#| (x = |A'(0+)|) is, however, 
super-universal, i.e. the same for all disordered elas- 
tic models. The significance of this cusp has been dis- 
cussed by several authors, ||, [l7| Notably, the fluc- 
tuations of the mean curvature of the pinning potential, 
-A£ =o (0) = (d%V{§)) 2 = +oo, which implies the ex- 
istence of multiple metastable minima in the effective 
potential. We will see further consequences below. In 
this T = regime, the neglect of iS' 3 ' appears consistent: 
solution of Eq. (fjj) gives <§4=o ~ ^(e 3 ) -C Rt=o> yielding 
only higher-order corrections. 

The convergence to this solution is non-uniform, how- 



ever, and for $ < eTi, the cusp singularity is rounded 
in a thermal boundary-layer. The boundary layer form 
follows from Eq. (0) for — through the ansatz 

R($) = -T?(e*r 2 r(0, with R"(0) = x 2 e 2 /(2_C - e), 
keeping cj> = e%$/Tj fixed of 0(1), x = %/e and Tj < e 2 . 
This gives r"(<f>) = 1 - v/l + 2 ,Jl5| matching the |0| 
cusp as <j) — > oo, but becoming smooth for <f> of 0(1). 

Even within the one loop Wilson RG, however, the 
neglect of the third cumulant is unjustified in the bound- 
ary layer. This follows since the 0[(R") 3 ] terms in 
Eq. (^) feed their boundary-layer scaling forms into . 
Naively balancing these terms with the rescaling part 
(~ — 2.S' 3 ') suggests ~ T l . This is, however, incon- 
sistent, since the remaining terms linear in would 

2 

then be O(Tj). Instead, the only consistent ansatz is 
5(3)($ 1; $ 2; $ 3 ) = ( ex )-2r 4 s (3)(^ l7 27 3 ). This in turn 
generates a feedback in Eq. (||) that is of the same order 
as the other surviving terms in the boundary layer. 

In fact, the scaling forms in Eqs. (|^-0) must be ex- 
tended to all cumulants. Defining U($) = AdAfU($Aj), 
then U($) = Ep,{ ai } h P <p ai ■ ■ ■ 4>a 2p + with 
/ap = A( 2 P)(0)/[(2p)!( X e) 2p ], where A( 2 f)(0) is the (2 P ) th 
cumulant of the random force (A( 2 )(0) = A(0)), and the 
boundary-layer replica vector <j>i = ex$/T/. This gives 
the non-trivial one-loop Wilson fixed point equation, 



P,i a i} 



Trin 



5 ab -d a d b u{$) =0, (8) 



3 



with scaling eigenvalues x^ — e — 2( of O(e) and X2 P 
■1 - 2p(0 - C) of 0(1) for p > 1. Expansion of Eq. (§) 
in number of replica sums gives coupled boundary layer 
equations for r and , and all higher cumulants. 

Unfortunately, the disappearance of e from Eq. (j|) 
(since X2 P f2 P ~ 0(1)), indicates that a one- loop anal- 
ysis may be insufficient. Worse, higher loop correc- 
tions to U appear naively extremely singular in 1/T: 
each additional contraction carries with it two 9$ deriva- 
tives, and hence two factors of 1/T, but only a sin- 
gle propagator proportional to T, so increasing powers 
of 1/T appear in higher loop diagrams. This makes 
the non-perturbative status of the boundary layer ex- 
tremely nontrivial. To verify its existence, we there- 
fore employ the Exact RG (ERG) method [|l8), taking 
from now on ( = for simplicity. The ERG extends 
the characteristic function £/($) to a non-local func- 
tional, W[$(r)], defined by a running effective action 
r(l) = (1/2)$ : Gr 1 : $ +Ui[$], with a cutoff Green's 
function (dUiq) = 5 ab (T/q 2 )(c(£) - c(^)). Then 

fkU0) = \ r &d l Gi : (G^-G^l+G, : ^T 1 ), (9) 



where 5's indicate functional differentiation. 



This 



equation can be expanded order by order in lii to 

generate a series expansion. The functional nature 

of U[<£>(r)] necessitates a further multilocal expan- 

sion = V; ..^r^-lv....!^:..;. 

in terms of translationally-invariant functions ("p-local 

terms") of p <&i and Xi variables. Remarkably, ex- 
it*) 

act integro-differential flow equations for the J7, can 
be obtained. The boundary layer ansatz, applied to 
each p-local term generalizes to l/j p ' ({$i}, {xi}) — 
e~ 2 h^ d Tiuf\{$i{\ , {A;Xi}), plus a nonlocal extension of 
the random- force terms (f2 P )- Inserting this into Eq. (|^), 
one finds that apart from the random-force terms, the 
left-hand-side is negligible for T; <C 1 and <fi ~ O(l). 
This gives a non-trivial hierarchy of integro-diffcrcntial 
equations (too complicated to present here) for the scal- 
ing functions u) with a well-defined limit as I — > oo. We 
have checked that the truncation to p < 2 is non-singular 
fL6| , i.e. the p = 2 term can be eliminated to obtain a 
single equation for the local function itj involving no 
divergent integrals. We believe this procedure, which ef- 
fectively resums the 1 /T divergences into a non-singular 
result, can be carried out systematically to all orders. 

This indicates the non-perturbative consistency of the 
boundary layer ansatz in the ERG. We now turn to some 
physical implications. A tree level calculation gives 



(($ 2 ) - ($) 2 )» ~ T 2n d\l a} U{Q) J d d r' G 2 "(r'), (10) 

where d 2r ^ yU is a combination of 2n th order derivatives 
of [/($), and G(r') = / d d q/q 2 ~ \r'\ 2 - d . Remarkably, 



inserting the boundary layer scaling form precisely re- 
produces the droplet behavior of thermal fluctuations, 
Eq. (||), with c„ ~ df" a} u{0)\ A further test of droplet 
predictions is that the simplest three-well quantity G3 = 
- $ h ) 2 ($ 6 - $ c ) 2 ($ a - $c) 2 } vanishes to 0{T). We 
have verified G 3 ~ 0(T 2 ) for d = 0, N = 1, 9 > 0. Some 
further results |l(| agree with recent calculations in this 
case where the simplest droplet picture is exact 

We now turn to the dynamics. Prior investigations 
]l5| have calculated the renormalization of the mean 
relaxation time at O(A), dfrj l = — A"(0)rj ; , leading 
through the boundary layer scaling A"(0) ~ 1/Tj to 
rj, = r;oexp[^(l/r ; - 1/T )] in Eq. (|). This is consis- 
tent with an energy barrier Ul — UoL e (with L — A ; _1 ). 

We begin our consideration of the distribution of 
relaxation times in a simplified calculation general- 
izing Eq. (^) to a random friction coefficient r\ — > 
rj(r), where the fluctuation-dissipation theorem requires 
(((r,t)((r' ,t')) = 2n(r)T5 d (r - r')S(t - t'). We take 
rj(r) with an independent and identical distribution P(n) 
at each r, with the characteristic function F(z) — 
— In J Q 00 drjP(i])e~ zv . In the standard Martin-Siggia- 
Rose (MSR) formalism, thermal and disorder averages 
over the equation of motion are calculated from a func- 
tional integral over $ rt , $ rt with the measure exp(— Sq — 
S st - Sd y ), where S = J rt -i$ rt V 2 $ rt , 

S st = Yt^T ■ ■ i$ rt JAW($ rtl , • • • , $ r J(ll) 

„ n -Jrt 1 -t n 



(12) 



and A (")(•) = (-1)"5 $1 •••9$„ {/(")(•). Setting = 
772: recovers the previous approach, but when F(z) is non- 
trivial Eqs. ( pryi2] ) comprise an unusual dynamical field 
theory. First consider the random friction model with 
A^™) = £A n ) = 0. Despite its highly non-quadratic ac- 
tion, this represents an FRG fixed manifold, with non- 
trivial response function moments, indexed by the "F- 
term" F(z), which remains un-renormalized. This fol- 
lows since Eq. (^]) becomes linear for zero pinning force 
/ = 0, hence the distribution of rj(r) is scale-independent. 

This changes drastically in the presence of pinning dis- 
order A( n ) ^ 0. Considering for the moment only the 
second cumulant n = 2, to O(A) F(z) is no longer scale- 
independent: diF t [z] = {dilnr] l )(zF{[z]+2z 2 F('[z]). This 
implies the rapid growth of higher connected cumulants 
of n, to wit F(z) =~J2q ^r-V iq) z q , with 

„ ^ e xM q )/T ie = ^) {%ho) ^) » W , (13) 

where a(q) = 2q 2 — g.Jlf], Thus all the renormahzed 
moments of the friction have the Arrhenius form, but 
since a(q) > q, they cannot be characterized by a single 
relaxation time ry ; . The quadratic form of a(q) obtained 



4 



in this approximation is consistent with a very broad, log- 
normal distribution of relaxation times r\. The more com- 
plete analysis below provides a consistent ansatz which 
extends this scaling to the full effective action, with, how- 
ever, non-trivial exponents a(g).pT[| 

To go beyond this simplified picture, we must con- 
sider all interaction vertices at a given "level" q involv- 
ing a fixed number q of total time derivatives (power 
of frequency lu 1 ), which includes e.g. F^)[0] but also 
an infinite hierarchy of higher cumulant functions, as 
in the static case q = 0. An examination of the one- 
loop low-temperature FRG shows that for each value of 
q one obtains a distinct hierarchy such that all terms 
grow proportionately to T^f 1 with a non-trivial exponent 
a (q). The hierarchies at different q are coupled in a spe- 
cial manner such that only levels q' < q feed into level 
q. This structure allows for independent exponents a(q) 
provided they satisfy the natural convexity constraint 
a(q) > a(q') + a(q - q') for 1 < q 1 < q - 1. Each 
member of the hierarchy scales as n\ q ^ times a distinct 
scaling function in the boundary layer regime <fi ~ 0(1). 
Coupled differential eigenvalue equations involving these 
functions determine the exponents a{q) by matching. 

We now illustrate this structure for q = 1, calculating 
the FRG for terms up to second cumulant order: 

Sdy = / rji^t^rt- iKtiiKt2$rtiG($rt-®rt 2 ) 
Jrt Jrt^ 



G 



rti*2*3 



where G(-$) = -G($). One finds that G contributes 
to the renormalization of 77: 



d l f}=(G (0)-A"(0))77, 



(15) 



?7 ; Af d A d 1 G($), correcting the naive 



where G($) - ^ 
A"(0) contribution to rj. Although G is nominally ir- 
relevant, diG = (—2 + e)G+ • • •, it supports a non-trivial 
boundary layer, in which G($) = ex<?(</>), satisfying 

= 2f"g + (l + f)g" - g'(0)(g - /') + 2/'(/"(0) + /") 

+3g'f + sg\(& <t>, 0) + 2/1100(0, 0, 0) + |/i O io(0, </>, 0), 

where /(</>) = -r"(0), H(^ 123 ) =r} l Af- 2d A^e 2 h(<l> 123 ). 
Matching to the well-behaved "outer" solution requires 
g{<j)) — >^->oo g+oo, a constant, while for cj) < 1, g(<j>) ~ 
g'(0)(t>- For given functions f,s^ 3 \h : this equation is an 
eigenvalue problem for g'(0) (indeed, analysis shows that 
g'(0) must be tuned to match the outer solution). From 
Eq. @, the exponent a(l) = g'(0) + /"(0) is thereby 
non-trivial as claimed. Moreover, a(l) is determined 
through the above equation in a manner depending upon 
the full boundary layer functions r, s, h. Similar consid- 
erations apply to q — 2 (details suppressed to save space), 
for which the one-time term is D J uu, and there are 3 



independent functions at the second cumulant level. The 
q = 1 functions feed into the q = 2 hierarchy in such a 
way that the a (2) exponent governing the growth of D 
is independent provided a(2) > 2a(l) self-consistently. 

In summary, we have described a consistent non- 
perturbative boundary layer structure in a static field 
theory for pinned elastic systems, successfully describing 
rare events dominating low-temperature thermal fluctu- 
ations. An extension to the dynamics leads to a broad 
distribution of relaxation times, for which the moments 
T n retain a (non-trivial) Arrhenius form. This condition 
would appear required for simple scaling arguments for 
creep to be even qualitatively correct, and is experimen- 
tally testable. Our investigations uncover an intimate 
connection between the highly non-perturbative bound- 
ary layer and droplet considerations, suggesting that an 
analytic approach encapsulating this physics might be 
successful in the near future. The more general applica- 
bility of these ideas to more complex glassy systems such 
as spin and structural glasses, as well as to the quantum 
regime, appears as a distinct possibility. 
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